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Abstract. We apply the generalized formalism and the techniques of the 
supersymmetric (susy) quantum mechanics to the cases where the super- 
\ potential is generated/defined by higher excited eigenstates (Robnik 1997, 

paper I). The generalization is technically almost straightforward but physi- 
cally quite nontrivial since it yields an infinity of new classes of susy-partner 
potentials, whose spectra are exactly identical except for the lowest n + 1 
states, if the superpotential is defined in terms of the (n + l)-st eigenfunc- 
tion, with n = reserved for the ground state. First we show that there 
are practically no possibilities for shape invariant potentials based on higher 
excited states. Then we calculate the isospectral partner potentials for the 
^ ■ following 1-dim potentials (after separation of variables where appropriate): 

(i) 3-dim (spherically symmetric) harmonic oscillator, (ii) 3-dim (isotropic) 
Kepler problem, (iii) Morse potential, (iv) Poschl- Teller type I potential, and 
(v) the 1-dim box potential. In all cases except in (v) we get new classes 
of solvable potentials. In (v) the partner potential to the box potential is a 
special case of Poschl- Teller type I potential. 
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1 Introduction 



In a recent paper (Robnik 1997, paper (I)) it has been shown that the for- 
malism of the supersymmetric (nonrelativistic) quantum mechanics can be 
applied also to the higher excited states (say, n-th state) of 1-dim potentials, 
generating new partner potentials isospectral to the original potential, ex- 
cept for the lowest n + 1 states which are simply just missing. There we gave 
the example of the 1-dim harmonic oscillator, which for all n > yields new 
classes of rational potentials. In this paper we present results of a straightfor- 
ward further application of this formalism to a few most important exactly 
solvable 1-dim potentials, namely (i) spherically symmetric 3-dim harmonic 
oscillator, (ii) 3-dim isotropic (spherically symmetric) Kepler potential, (iii) 
Morse potential, (iv) Poschl-Teller type I potential, and (v) 1-dim box po- 
tential. 

Following the seminal papers of Witten (1981) and Gendenshtein (1983) the 
methods of supersymmetric (susy) (nonrelativistic) quantum mechanics have 
quickly developed and it has been realized, that (1) there exist partner po- 
tentials with precisely the same energy spectra except for the ground state 
(n = 0) (whose wavefunction 0(x) = ipo(x) is used to generate/define the su- 
perpotential W(x) - see below) 0, and that (2) if they are "shape invariant", 
their spectra and wavefunctions can be exactly and analytically solved. It is 
believed that the list of such susy-0 shape invariant partner potentials is now 
complete and finite (Levai 1989, Barclay et al 1993), and therefore quite lim- 
ited in use. The research has been later further developed also in direction 
of applying the WKB methods to such classes of Hamiltonians, including 
the search for improved simple quantization conditions which would be ex- 
act in case of susy shape invariant potentials (Barclay, Khare and Sukhatme 
1993, Barclay and Maxwell 1991, Barclay 1993, Inomata, Junker and Su- 
parmi 1993, Junker 1995, Robnik and Salasnich 1997), and also in direction 
of exploring the applicability of the path integral techniques (Inomata and 
Junker 1991,1994). One of the nicest presentations of susy quantum mechan- 
ics was published by Dutt, Khare and Sukhatme (1988), henceforth referred 
to as DKS. We will use their notations. It should be mentioned at this place 

3 the ground state energy Eq is missing in the partner Hamiltonian H + , so that its 
groundstate Eq — E^ 
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that the ideas involved behind the susy property and shape invariance were 
formulated first by Infeld and Hull (1951), where they were called the "fac- 
torization method", and these authors refer further to the related ideas in 
the works of Schrodinger (1940,1941). 

Thus we shall use the notations of DKS, employed also in (I), and present the 
brief outline of the susy-m formalism in the section 2. Then we demonstrate 
in section 3 that there are practically no cases of susy-m shape invariance for 
n > 0. In sections 4-8 we present the partner potentails of (i) through (v), 
and in section 9 we draw the general conclusions and discuss the results. 



2 Generalized supersymmetric formalism 

The main point of this section is to briefly review the formalism of the susy- 
n quantum mechanics, which thus can be generalized to arbitrary higher 
excited eigenstates <f>(x) = ip n (x), n — 0,1,2,..., used to generate the 
superpotential W(x), namely 



where (f>'(x) = d<f>/ dx, \i is the mass of the particle moving in the V~ potential, 
lnh, is the Planck constant and n is the quantum number equal to the number 
of nodes of the eigenfunctions ip n (x) of the starting potential V~(x). The 
energy scale is adjusted so that the (n + l)-st energy eigenvalue is exactly 
zero, E~ = 0. The corresponding Hamiltonian is H = — + V~(x), 
and the Schrodinger equation reads 

H~r n = H-<t> = {~yJ^ 2 + V-{x))4> = 0. (2) 

Obviously, because 4>'{x) ^ at the nodes yj, the superpotential W(x) will 
have singularities at the nodes yj, j = 1, 2, . . . , n of (f>. However, this does not 
invalidate our derivation, but it merely means, as will become clear later on, 
that the partner potential generated by diverges to +oo when x — > yj, for 
any j = 1,2, ... ,n. This implies that the potential wells are well defined be- 
tween two consecutive singularities and that they do not communicate with 
solutions in the neighbouring wells. Thus if n = we have the common case 



4 



of usual susy potentials defined on (—00, +00), if n = 1 we have two sepa- 
rated potential wells, each of them on a semiinfinite domain, for n = 2 we 
have one infinite potential well on a finite domain between two nodes y\ and 
U2, and two binding potential wells on the two semiinfinite domains (—00, y\] 
and [y 2 , +00), and so on. The (partner) potentials constructed in this way 
are nontrivial and certainly very interesting since they contribute to our list 
of solvable potentials which now becomes truly very rich and infinite in its 
contents. 

In order to make this paper selfcontained we will build up the formalism nec- 
essary to construct the partner potentials and to define the shape invariance, 
following DKS, in order to demonstrate that the susy formalism does not 
break down anywhere on its domain of definition, and to define the language 
needed to talk about further results that we shall present in this contribution. 

First we express the starting potential V~(x) in terms of the (n + l)-st 
eigenfunction <p(x) = ip n {x), by solving (0) 

which is regular everywhere, because at the nodes yj the second derivative 
4>"(x) = cPip/dx 2 also vanishes with 0. Thus the basic Hamiltonian H~ reads 

h 2 d 2 6" 

H - = h { -l> + J ] - (4) 

The two important operators are: 

At = 4=(-4--l), (5) 

and 



2fi dx 



which gives 



/2/7 dx 

H~ = A^A. (7) 
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We further define the partner Hamiltonian H + and the partner potential V + 
as 



where 



or 



V+( X ) = _V-( X ) + ^!(^)2. (10) 

The potentials V + and V~ are called susy-n partner potentials. We will 
show that they have the same energy levels, except for the (n + 1) lowest 
states of V~ for which there are no corresponding states of V + , so that the 
ground state of the latter one is Eq = En+i- All higher states have then 
identical energies. From equation ( |T0"D we see explicitly that at every node 
yj, j = 1, 2, . . . , n of the defining eigenstate = ip~ the partner potential(s) 
will have a singularity of the type l/(x — yj) 2 such that V + (x) — > +oo when 
x — > yj, so that every branch of the partner potential will be confining up to 
infinity, and the solutions in various branches do not communicate. Thus for 
each n we shall find (n + 1) (branches of the) partner potentials. 

In terms of the superpotential W defined in equation (|1|) we can write 



(f)(x) = = exp (-^^ / W{x)dx), (11) 



h 

which is well defined in the definition domain of any of the branches of the 
partner potential, and obviously will go to zero where W has the poles 
l/(x — yj), as it should happen. 

Some of the relationships can be rewritten/reformulated in terms of the su- 
perpotential W(x) now: 
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(12) 



Further we observe, 



V ± (x) = W 2 (x) ± JL;W'(x), W\x) = (13) 



and also 



The commutator of the operators A and is 



1 ' J dx v ; 



Now we have all tools at hand to show that the susy partner potentials V~ 
and V + are isospectral except for the lowest (n + 1) states of V~ which have 
no counterpart in V + , so that its ground state is Eq = E~ +1 . 

The demonstration, following DKS, is very easy: First we find that if ■?/>" 
is an eigenfunction of H~ with the eigenenergy then Aip^ is an eigen- 
function of H + with the same energy: 

H + (M~) = AdAtl,- = AH~r m = AE~ m r m = E~A^~. (16) 

Now we show that this applies only to the eigenstates m higher than n, m = 
n + 1, n + 2, . . ., by considering the normalization condition, by writing the 
normalized state = C m Aip m , and calculating the normalizing coefficient 

Cm, 



||^+|| 2 = C 2 m < A1>~\A1>- >= C 2 m < ^\A + A^ m >= C 2 m E-U-\\ 2 . (17) 

If all ■?/>" are normalized (they are certainly orthogonal, because we deal with 
one dimensional systems, where degeneracies are forbidden due to the Sturm- 
Liouville theorem (Courant and Hilbert 1968) and therefore all eigenstates 
must be orthogonal), then 
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which implies that the construction succeeds only iff > 0, implying that 
m > n. Thus the two Hamiltonians H~ and H + defined in (§) and in (||) 
are isospectral except for the lowest (n+ 1) eigenstates of which have no 
counterpart in H + . 

Counting now the eigenstates of H + from m = 0, 1, 2, . . ., where m = is 
the ground state, and m is the number of nodes of the (now also normalized) 
eigenfunction we have 

Of course it is easy to show that, conversely, for every eigenstate ip m °f H + 
there exists the normalized eigenstate of H~, namely 

^- +l+m = ^=A^^ m = 0,1,2,... (20) 



This completes our proof of isospectrality, generalized to the case that the 
generating function (ft of the superpotential W, defined in equation (|l|), is a 
higher excited wavefunction, namely <fi = n = 0, 1,2, . . .. As we have 
seen, the formalism of superpotential and of the partner potentials works 
everywhere except at the singularities located at the nodal points yi of 4>, 
where the partner potential V + goes to infinity as l/(x — yi) 2 , thereby defin- 
ing several branches of V + well defined on their disjoint domains of definition. 

We have demonstrated that if one of the partner systems (the Hamiltonians) 
can be solved completely (by calculating the energy levels and the eigen- 
functions), then the susy formalism enables one to solve the partner problem 
completely, following equation flT5|). One of the most important cases is of 
course the harmonic oscillator, which has been discussed in (I), whilst in this 
paper we study five more examples (i)-(v) announced in Introduction. 
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3 Is there shape invariance for n > 0? 



If the solutions for the two partner Hamiltonians are both unknown, then 
another approach is necessary to solve them. In case of the standard susy 
formalism with n = we have the important class of the shape invariant 
potentials. As is well known (DKS) the shape invariance of the two partner 
potentials V~ and V + is defined by 

V + (x;a ) =V-(x;a 1 ) + R(a 1 ), (21) 

where ao is a set of parameters, a\ = /(ao) and R(ai) is independent of x. 
The procedure is now (essentially embodied in the factorization method of 
Infeld and Hull (1951)) the following. Consider a series of Hamiltonians H^ s \ 
s = 0,l,2,..., where H <°) = H and if W = H + , by definition 

H {s) = + V-(x; a s ) + £ R(a k ), (22) 

where 

a a = f'(a ) = fo...of(a ). (23) 

s 

Now compare the spectra of if ^ with H^ s+1 \ and find 

= ~^-^ + V ~ {x;as+l) + ^ R{ak) > 
H {s+1) = -f^- 2 +V + (x;a s ) + ±R(a k ), (24) 

* k=l 

Thus it is obvious that and if ( s+1 ) are susy partner Hamiltonians, and 
they have the same spectra from the first level upwards except for the ground 
state of whose energy is 

E^ = J2R(a k ). (25) 
k=i 

When going back from s to (s — 1) we reach = H + and = H , 
whose ground state energy is zero and its m-th energy level being coincident 
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with the ground state of the Hamiltonian H^ m \ m = 1, 2, . . .. Therefore the 
complete spectrum of H~ is 

rn 

E- = J2 R M E =0. (26) 
fc=i 



The generalization of shape invariance to the case of any n > is straight- 
forward, but it results in higher complexity and therefore it is more rarely 
satisfied by the specific systems. By repeating the above argumentation we 
reach the conclusion that, when fl2"T|) is satisfied for a superpotential W with 
given n, then we cannot calculate the entire spectrum of the shape invari- 
ant potential/Hamiltonian H~ , but only the subset (subsequence) of period 
n + 1, namely 

m 

K +m{ n + i) = E RM, E~ = 0, m = 1, 2, . . . . (27) 

k=l 

In the special case n = we of course recover the formula (p6|) . For n > we 
have none example of susy-n shape invariance so far. In fact we shall show 
now that there is no susy-n shape invariance for n > 0, unless we find some 
rare exceptions. 

In fact we can see in equation ([H]) that if V~(x;ai) has no singularities 
as a function of x, then the partner potential V + (x;ao) also has no singu- 
larities, because R{a\) is just a constant and independent of x. Therefore, 
for n>0a potential V(x) cannot be shape invariant if it has no singulari- 
ties. The only possibility then is that V~(x; a x ) has singularities at the same 
places yi as the partner potential V + (x; a ): This, however seems also impos- 
sible, because for n > the partner potential V + obtains new singularities 
between those of V~ , corresponding to the nodes of the ra-th eigenfunction 
of V~ . Therefore, unless we find some pathological exceptions, it seems that 
there is no shape invariance for higher excited states n > 0, at least not in 



the sense of the definition (|2l|). Perhaps some other functional relationships 



yet to be discovered might lead to some other type of "shape invariance". 
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4 The 3-dim spherically symmetric harmonic 
oscillator 



Let us consider a few examples of susy-m partner potentials, first the 3- 
dim spherically symmetric harmonic oscillator. (The case of 1-dim harmonic 
oscillator was calculated and discussed in (I).) To prepare some generali- 
ties of 3-dim spherically symmetric potentials V(r) we first write down the 
Schrodinger equation Hip = Eip, namely 

^& + ll> + ^ + v ^'- E ^ < 28 > 

where /i is the mass of the particle and H, L 2 and L z are the usual notations 
for the Hamilton operator, (square of the) angular momentum, and the z- 
component of the angular momentum. The quantum numbers of the latter 
two will be denoted by I and m, so that / = 0, 1, 2, . . . and m — 1,1 — 1,1 — 
2, . . . , —(I — 1), —I, and the eigenvalues of L 2 are /(/ + l)h . Due to the 
spherical symmetry we have thus the separation of variables 

il>{r,6,<p) = R(r)Y lm {6,<p), (29) 
where YJ m are the spherical harmonics, so that 



h 2 . d 2 2d, . > 1(1 + l)h 2 , 
W^- + rTr ) + V{r) + ^- ]R - ER ' (3 °» 
and after the substitution 



R{r)= X {r)/r, V t (r) = V(r) + 1(1 ± ^ - , (31) 
the radial Schrodinger equation becomes finally 

r h 2 d 2 



l , r > rV l (r)} X (r) = E X (r), (32) 

which is now just a 1-dim problem with the effective potential Vi(r) defined 
in equation (|3l|) , for the physical range of definition r > 0. From equation 
29T) follows the normalization condition 
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oo 

2/ 



dr X (r) = 1. (33) 
Now we look at the specific case of the 3-dim harmonic oscillator defined 

by 

V{r) = -/MJ 2 r 2 . (34) 

Introducing the constants 



, _ 2fiE 2 fiu 



e = ^, A — (35) 



we rewrite the Schrodinger equation ( |32]) as 



(// ., + ^ - - = 0, (36) 

with the solutions 

Xn (r) = Cr l+1 e^ x2r2 F(-n, I + ?; AV), (37) 

where n is the radial quantum number n — 0, 1, 2, . . . (in denoting the eigen- 
functions and eigenvalues we shall suppress the quantum number I for brevity 
of notation), and the energy eigenvalues read 

3 

E n = (2n + Z + (38) 

In the above -F(a, b; x) is the confluent hypergeometric series, [] and the nor- 
malization constant C in (|3T| ) is equal to (Goldhammer 1963) 

2;+ 2-n (2 ^ + 2n+1)! , 

[ 0Fn![(2/ + l)!!] 2 J A ' ldyj 

The superpotential defined according to (fj) using the n-th eigenfunction ([37]) 
reads 



4 It is also called Rummer's function and denoted by M(a, b, z) e.g. in Abramowitz and 
Stegun (1965) p.504. 



12 



" l J v^ 1 r L + (/ + 3/2)F(-n,Z + 3/2;AV 2 ) JI 1 J 

so that for the susy-0 (ground state) superpotential we have 

W (r) = JEur - {l+ r l )H . (41) 

Now we can calculate the partner potential starting from the shifted 
potential 



1 2 2 l(l + l)h 2 

-au 2 r z + — 

2 P 2fir 

where we see E~ = 0, and thus obtain the partner potential 



V-(r) = -imJ 1 ,' 1 + X Z.T ~ ( 2n + ' + 3 / 2 > fiw ( 42 > 



2/i <ir 

1 2 2 (/ + l)(/ + 2)^, 2 /rt , V 

= +- (2n + Z + -)nw + / n , 

2n^ 2/io; 2 (1-n) n ^ 2 (l) n 

(Z + S^i^O) 1 1 ; h [ (/ + 5/2) 1 ; (Z + 3/2) F(0) Ji 

(43) 

where we used the short-hand notation for the confluent hypergeometric func- 
tion 

F(i) = F(-n + i,l + ^ + i-\ 2 r 2 ), i = 0,1,2. (44) 

It is readily seen that for n = we get the shape invariant case of the two 
partner potentials 



(r, I) = r + 2/ir2 ~ ( Z + ^ 
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T,+ , A 1 2 2 (l + l)(l + 2)k 2 n 1.. 

Wr,i) = + - ^ r2 + ^ 

K+(r,/) = ^T(r,i + l) + 2^o; 

a = I, a\ — I + 1, -R(ai) = 27ia>. (45) 

As the last point of this section we give the explicit expressions for the two 
lowest excited states (n = 1,2): 

N 1 2 2 (/ + l)(/ + 2)/i 2 5_ 

= 2^ V + - ^ r2 -(/ + -)^ + /i 

Ji = 2/? ^(l -z)- 1 [l + 2^(l-^)" 1 ] 



(46) 



and 



I ,-(/■) = 2^ + - 2V -fl + ^ + fr 



x {1 - fa + 2z[-/3x + 2/3 (l - 2/3 z + fafcz 1 )- 1 ^ - ft*)]} 

(47) 



The wavefunctions can be calculated using formula ([19|), where we have 

A = ^s + Hr -M- (48) 



and using the eigenfunctions fl3T|) and ( |4"0| ) we can write down the expression 
for the eigenfunctions of the susy-n partner potential, namely (bearing in 
mind x ~ = Xn, horn (§7|)), 

Xp = I =AXn+X+p> (49) 

^n+X+p 
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with (using the equation (|38])) 

K+i+p = E n+i+p -E n = 2(p + l)hw, p = 0, 1, 2, 



(50) 



5 The 3-dim spherically symmetric Kepler prob- 
lem 

The 3-dim spherically symmetric Kepler potential is 



V(r) = —. 

r 

When using this in equations ( ^ - |3^) we find the solutions 



Xn{r) 



Cr l+1 e-? r F(-n,2l + 2;\r), 



(51) 



(52) 



where n = 0,l,2,...is the radial quantum number (counting the number of 
nodes of Xn( r )) an d the energy eigenvalues are 



E„ 



2h 2 {n + l + l) 2 ' 
with the definitions of the parameters 



A 



(n + l + l)a' 
The normalization condition is 



fie 2 



X 2 n (r)dr = 1, 



so that the constant C is equal to 



C 



(N + l)\ 



j5/2jV3(2Z + l)!\| (N-l-l)V 
The susy-n superpotential as defined by (jl]) reads 

h l + l A nAF(l) 



N = n + I + 1. 



WJr) 



2fi L r 2 (2/ + 2)F(0) J 



(53) 



(54) 



(55) 



(56) 



(57) 
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whose special susy-0 case is 

w " (r) = V 2 (TTTjA " < 58) 

In the above we used a similar short-hand notation for the confluent hyper- 
geometric series (function) as in the previous section, namely 

F(i) = F(-n + i,2l + 2 + i;\r), i = 0,1,2. (59) 
Further we calculate the partner potentials 



. . . , „ e 2 Z(Z + l)ft 2 /ze 4 

r + 2^ + 2fi 2 (n + ' + l) 2 
nAW .fl-n) nFj) 



/i(2/ + 2)F(0) L (2/ + 3) v 7 (2/ + 2)F(0) J 

Now in case of susy-0 partner potnatials we recover the shape invariance 
property embodied in the following relationships 



V-(r) ~ e2 | ^ + 1 ^ | ^ -V-(rl) 
V+(r) - (/ + !)(< + 2)fi 2 ^ 

V + (r,l) = V -(r,l + l) + ^[- 1 ' 1 



2/i 2L (Z + l) 2 (/ + 2) 2J ' 
V + (r,a ) = y _ (r,ai) + i?(ai), a a = l, a x = I + 1, 
^ x /ue 4 r 1 1 , 



2a 2L (/ + l) 2 (/ + 2) 2r 

(61) 

Finally we calculate the partner potentials for the two lowest excited states. 
We find 
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v+(r) _ e 2 (l + m + 2)g ^ 

r + 2^r 2 + 2h 2 (l + 2) 2 ^ h ' 
\ 2 h 2 r . Ar 



-2 



fi(2l + 2) 2L (2l + 2y 



(62) 



and 



V2[T) ~ 7 + 2^ 2 + 2h 2 (l + 3) 2+h ' 
8/ie 4 /3 



r (l -2/?0z + (3 (3 1 z 2 )' 1 x 



{i + 3) 2 h 2 

x [-ft + z(3 (l - f3 lZ ) 2 (l - 2(3 z + Pafcz 2 )- 1 }, 

2r 2fie A r 



(3 l = 21 + 2, /? 1 - 1 = /5- 1 + l, 



Finally we compute the eigenenergies 



(l + 3)a {l + 3)h 2 ' 

(63) 



4 

E; = E~ +1+p = E n+1+P -E n = ^[(n+l+l)- 2 -(n+p+l+2)-% E+ = E~ +1 . 

(64) 

and the eigenf unctions using equations (ffl), (fl9|) , (|52|). 



X}{r) = 7= Ax ra+1+P (r), X„+i+ P (^) = Xn+i+ P (r), p = 0, 1, 2, . . . 
'E n+ i +p 

(65) 



6 The 1-dim Morse potential 

In this section we analyze the supersymmetric aspects of the Morse potential, 
defined as 
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V(x) = A(e- 2ax -2e~ ax ), (66) 

which is an important model potential. It has a finite number of eigenstates. 
The notation for the constant A here should not be confused with the (ladder) 
operator A from section 2. The Schrodinger equation 

*^ + ^(E-V(x)Mx) = 0, (67) 

can be rewritten as 

d 2 i> I dip , 1 n + s + \ s 2 x n 

where we use the notations 



Now we try the Ansatz 

such that the function u(£) must satisfy the simpler equation 

£u" + (2s + 1 - f )u' + nu = 0, (71) 
having the solution in terms of the confluent hypergeometric function 

u(0 = F(-n,2s + l;0, (72) 
so that finally the explicit solution obtains the form 

MO=CCe^ /2 F(-n,2s + l;0, (73) 
with the energy spectrum 



* = -^-^<» + 5>l* » = 0,1,2....,<(^-1) (74, 
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and the quantum number n runs up to the maximum value. 

Now employing the same procedure as before we find the susy-n super- 
potential W n (x), namely 

„, , , h r 1 . na£F(l) . 

W'^^H ^ 1 (75) 



whose special and well known ground state (n = 0) value is 

W (x) = + | - V^e" 08 ". (76) 

Again, we use here the short-hand notation for the confluent hypergeometric 
function 

F(i) = F(-n + i,2s + l+i;£), £ = 0, 1, 2. (77) 
The starting shifted potential of flBTSD is 



V~(x) = V(x) -E n = A(e~ 2ax - 2e~ ax ) + A[l - -^=(n + -)] 2 , (78) 



V2/Z4 V 2' 



and we get the isospectral partner potential 



V+(x) = Ae- 2ax -(2A- x ^ah)e- ax + A[l-^=(n + l)] 2 + I n 

V [i y/2fj,A 2 

In = [F( i) i (IzlOgg) i < F2 W i 



(2s + 1)//F(0) L v ' (2s + 2) (2s + l)F(0) J 

(79) 

The shape invariance is recovered for n — 0, namely the shifted starting 
potential is 



V -(x) = V(x)-E = A(e- 2ax -2e-n+A(l-^==) 2 = V ~(x; 2A), (80) 
and the isospectral partner potential 
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V+(x) = Ae- 2ax -{2A-^ah)e- ax + A{l- 1 ^=f = V+(x;2A), (81) 

where the notation Vq(x;2A) implies dependence on the coefficient 2 A in 
the Morse potential (|66|) . Then we find the shape invariance relationship 



Vff(x;a ) = V (x; ax) + R(ax) 
a = 2A, 



+ 

a\ = 2A — \ — an, 

V ^ 



2^/2A^' 



(82) 



The isospectral partner potentials of V 1 and of V 2 are calculated in a 
straightforward manner: 



V+(x) 
h 



<2A 



3ah 



Ae~ 2ax - (2 A - \l —ah)e~ ax + All 

V V H ' 2^A' 

a 2 f3 hl; 



+ h, 



;i-/3oO _1 [i + M(i-AO _1 ] 



(83) 



and 



V+(x) = Ae~ 2ax - (2A - J —ah)e- ax + A(l - —j==) 2 + I 2 

y /i 2^2/j,A 

h = 2a 2 f3 h 2 ai-2M+PoPier 1 x 
x [i - fat + 2/3 e(i - 2M + PoPier 1 ^ - m 2 \ 



where we use the notation for the parameters Po, Pi, 

A) = 2s + 1, A" 1 = /5o _1 + 1- 



(84) 



(85) 
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Finally the normalized wavefunctions of the partner potential V^{x) are 
calculated easily 



Ipp — J — = Alp n+l+p , V'n+l+p ~~ i^n+l+pi E n+ i +p — E n+ i +p E n , 

E n +i-\- p 

(86) 

where again p = 0,1,2,..., and A here is operator defined in (0) with the 



superpotential given in (|75|). 



7 The 1-dim Poschl- Teller type I potential 

The 1-dim Poschl- Teller potential is defined as (Fliigge 1971) 



l„ t A(A-l) , B(B-1 

t; v o[ — r~o 1 5 — 

2 sin ax cos^ ax 

where 



W= H[^ + S (87) 



h 2 a 2 n , , 

V = , < ax < -. (88) 

The solution is well known 

4>n(x) = C sin A ax cos B axF(—n, A + B + n; sin 2 ax), (89) 
with the energy spectrum 

E n = ^V (A + B + 2n) 2 , n = 0,1,2,.... (90) 
The susy-n superpotential is 



ah 2n(A + B + n)smaxcosaxF n (l) 
W n (x) = == \A cot ax — B tan ax — m=n m 

afi 

Wq(x) = t=(A cot ax — B tan ax) 

v 2/U 



(91) 
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where we employ the short-hand notation for the hypergeometric series (Abramowitz 
and Stegun 1965 p. 555) 

F n (i) = F(-n + i, A + B + n + i, A + - + i; sin 2 ax). (92) 

2 

The starting shifted potential is 

V~(x) = V(x) - E n = \v [ A{A ~ 1] + B{B ~ l) -(A + B + 2n)% (93) 

2 sin ax cos 2 ax 

and the isospectral partner potential 

2 sin ax cos 2 ax 

2V n(A + B + n) r „ , w 9 . 

= (V + 1)F.(0) {F - (1)(1 - 2m " M!) + 

2(-„. + 1)(A + g + n + 1)F.(2) 2n( 4 + g + ■ „)fg(l) 

1 (A + 3/2) (A + 1/2)F„(0) 1 ' 

(94) 

In case n = we again recover the shape invariance property 

V -(x) = \v [^^ + ^f^-{A + BY\ = V Q -{x,A,B), 
2 sm ax cos 1 ax 

V + (x) = l Vo [ A (^ + B l^-(A + Bn = V+(x,A,B), 
2 sm ax cos z ax 

V + (x,{a }) = V -(x,{ ai }) + R({ ai }), 

{a } = {A, B}, { ai } = {A + l,B + l}, R({ ai }) = V (A + B + 2). 

(95) 

The isospectral partner potentials for the lowest two excited states n = 
1, 2 are 

vra = V ol d<4±i) + £(£±i)_ ( . 4 + B + 2)2| + /l , 

2 sin ax cos 2 ax 
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2V (A + B + 1) 



x 




} 



(96) 



and 



V+(x) 




(A + B + 4f] + I 2 , 



+ [ 



2(A + B + 3) 4(A + B + 2)F$(l) 
(^4 + 3/2) + (A + I)F 2 (0) 



] sin 2 ai cos 2 ax} 



(97) 



The procedure to obtain the normalized eigenfunctions of the isospectral 
partner potential is the same as in the previous sections and in general 
it is explained in the section 2. 

8 The 1-dim box potential 

The 1-dim box potential is defined as 

V{x) = for - L/2 < x < +L/2, and V(x) = oo otherwise. (98) 
The solutions of the Schrodinger equation 



2fidx 2 



+ (E — V(x))ip = 0, 



(99) 



are 





cos 



nirx 



L 



odd n — 1, 3 
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, . . 2 . nwx 

Vn\x) = y^ sm— even n = 2,A,... 

(100) 

with the energy spectrum 

h 2 ir 2 9 

n= 2^ ' 7l = 1 ' 2 .---- ( 101 ) 
The shifted starting potential is 

V^~(x) = — _E n for — L/2 < x < +L/2, and V"(a;) = oo otherwise, 

(102) 

the superpotential is 

W n (x) = — j|=^, (103) 



and therefore using equation (JTJ|) we get the partner potential 



Ti 2 1 

V+{x) = -E n + — - , odd, 

H cos z k n x 

h 2 1 

V+(x) = -E n + — 



/i sin k n x 
L ' En= 2~il' 



mi h 2 2 

fen — "i - j -C^n = ~^~k , n = 1, 2, o, 



(104) 



Thus the isospectral partner potential to the 1-dim box potential is a special 
case of the Poschl- Teller type I potential, defined in (fffi) . The eigenenergies 
for V+(x) are 



= E n+l+p -E n = Yjz[{n + l+p) 2 -n 2 }, p = 0, 1, 2, . . . (105) 
and the eigenf unctions are easily obtained by applying the formula (IT 
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9 Discussion and conclusions 

In this paper we have applied the supersymmetric formalism of the (nonrel- 
ativistic) quantum mechanics, introduced and explained in (Robnik 1997), 
to the higher excited states of five different specific exactly solvable poten- 
tials, namely (i) 3-dim (spherically symmetric) harmonic oscillator, (ii) 3-dim 
(isotropic) Kepler problem, (iii) Morse potential, (iv) Poschl- Teller type I po- 
tential, and (v) the 1-dim box potential. In all cases except in (v) we get new 
classes of isospectral partner potentials which thus also fall into the class of 
exactly solvable potentials although they can be and typically are quite com- 
plex. The most important case of the 1-dim harmonic oscillator was treated 
in detail in paper (I) (Robnik 1997). 

We have also shown, in section (3), that for higher excited states there is 
generally no shape invariance of the usual type, and it remains to be investi- 
gated if there are some other types of "shape invariance", generated by some 
other functional relationships between the starting and partner potentials, 
such that they would allow for an exact solution of the problem. At present 
we do not know any specific cases of susy-n shape invariance with n > 0. 
Further calculations for susy-n partner potentials isospectral to some other 
well known exactly solvable potentials, not analyzed in this paper and in (I), 
remain as a future project. 
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